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Abstract 

We consider the transitions from the Z{,(10610) and Zf,(10650) resonances to the 
states of bottomonium with emission of a pion. The Zb resonances are viewed as 
'molecular' objects of a large spatial size made of heavy mesons, while the states 
of bottomonium are considered to be compact, so that an application of the leading 
(dipole) term in the QCD multipole expansion is assumed to be justified. In this way 
we calculate the ratios of the decay rates to the final states T(nS) it with n = 1,2, 3, 
and the ratio of the decay rates to hb(kP) ir with k = 1 and 2. We find our estimates 
in a reasonable agreement with recent experimental data. 



The isovector 'twin' resonances Z^ = Z&(10610) and Z' b = Z;,(10650), recently found 
by the Belle Collaboration pQ near the respective thresholds B*B and B*B*, are naturally 
interpreted (2j |3l HI [5] as molecular states made of the heavy meson-antimeson pairs. Each 
of the two new resonances is observed through the decay to either T(nS) tt with n = 1, 2, or 
3, or to hb(kP) tt with k = 1 or 2. Moreover, the decays to the states of ortho- and para- 
bottomonium are found to have comparable strength with no suppression of either of them 
by the heavy quark spin symmetry. This behavior is natural within the interpretation of 
the Zb resonances as being molecular S'-wave states of the heavy mesons: Zb ~ B*B — B*B 
and Z' b ~ B*B*, since the total spin of the bb quark pair within a meson system is not 
fixed |2j. Although at present the type of the threshold singularity in the heavy meson - 
antimeson channel (bound, virtual, or resonant state) corresponding to the Z b peaks is not 
known, it appears clear that the Z b peaks result from a strong dynamics of very slowly 
moving mesons near the threshold. The notion that the heavy b and b quarks at the Zb 
resonances are moving at distances longer than the characteristic size of bottomonium is 
also in a qualitative agreement with the recently available data [6] on the yield of different 
radial excitations of bottomonium in the decays of these resonances to T(nS) tt and hb(kP) tt. 
Namely, the yield does not diminish with the number of excitation in spite of kinematical 
suppression for production of heavier states. This implies, at a qualitative level, that the 
overlap of the bottomonium states with a widely separated heavy quark pair in the initial 
state increases with the excitation number due to larger spatial size of the excited states. 

The purpose of this paper is to suggest an approach to calculating the relative rates of 
decay of the Z b resonances to various radial excitations of bottomonium with emission of light 
mesons, and thus to quantify the theoretical estimates of the relative strength of the observed 
transitions Z\ — > T(nS)n and Z b — > hb(kP) tt. We assume that the bottomonium bb 
system in the final state is pure color singlet and is sufficiently compact, so that its interaction 
with soft gluon field can be considered within the multipole expansion in QCD [7J[8] with the 
leading term being the chromo-electric dipole. The transition of the heavy bb pair from the 
initial 'molecular' state to bottomonium is due to this interaction at short distances whose 
scale is set by the bottomonium size, while the (soft) gluon field induces the transition of 
the light quark- ant iquark-gluon components of the initial 'molecule' to the light hadron(s) 
in the final state. In this picture the specific form of the heavy bb 'overlap' amplitude is set 
by the wave function of the initial state, the chromo-electric dipole interaction, and the wave 
function of the bottomonium state. Therefore, given a model for the latter wave function 
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for various radial excitations, one can evaluate the relative strength of the transitions to 
those excitated states of bottomonium. In our estimates in this paper we use the simple 
model with the Cornell potential [9]. As for the initial state wave function of the bb pair we 
use the short- distance part of that for a slowly moving pair. Moreover, the chromo-electric 
interaction links the color singlet finite state to a color-octet initial bb pair. Clearly, such 
state is present in a molecular heavy meson-antimeson system. Indeed, in the colorless B or 
B* mesons the color of the b quark is correlated with the color of the light antiquark. Then 
in a well separated meson-antimeson system the color of b in the meson is fully uncorrelated 
with that of b in the antimeson, so that a color-octet bb pair is present with the statistical 
weight 8/9 and the statistical weight of a colorless bb state is 1/9. At short distances there 
is a weak (Coulomb-like) repulsion between b and b in the color octet state. Although 
suppressed by the color factor l/(N^ — 1) = 1/8 the effect of this repulsion is noticeable 
at a small momentum of the heavy quarks, and we take it into account. Furthermore, it is 
important for the discussed approach that, even though the orbital angular momentum of 
the heavy mesons in a molecular system can be fixed (5* wave in the Z b resonances), the 
orbital angular momentum of the bb pair in such system is generally not fixed (2j [3] due to 
the motion inside the mesons. Thus the chromo-electric dipole transitions to the S"-wave 
T(nS) of bottomonium proceed from the P-wave state of the initial bb pair, while those 
transitions to the P-wave hb(kP) levels are dominantly from the initial bb S-w&ve pair, since 
the wave function for a P-wave is suppressed at short distances. 

In the calculation in this paper we use the Hamiltonian for the chromo-electric dipole 
interaction in the form 

H E1 = -^Cf-E a (Q) , (1) 

where £ a = — ^2 ^ s the difference of the color generators acting on the quark and antiquark 
(e.g. tf = X a /2 with A a being the Gell-Mann matrices), r is the vector for relative position 
of the heavy quark and the antiquark. Finally, E is the chromo-electric gluon field strength. 
We therefore write the amplitudes of the discussed decays in the form 

(T(nS) ir\H E1 \Z b ) = C s A nS E n (Z ■ f) 

(h b (kP)n\H El \Z b ) = C P A kP (p w ■ [Z x h}) , (2) 

where E v (p n ) is the pion energy (momentum), Z, T, and h are the polarization amplitudes 
of the initial and final resonances, and Cs, Cp are constants that do not depend on the 
excitation number of the final state, while this dependence is contained in the amplitudes 
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A n s and A k p describing the overlap integrals with the dipole interaction ([Q): 

A nS = ! R nS {r)rRf{r)r 2 dr , (3) 



A kP = J R kP (r)rRf\r)r 2 dr . (4) 

In the latter expressions R n s {Rkp) are the radial wave functions of the bottomonium S (P) 
wave states and Rg (Rp ) are the radial wave functions of the color-octet bb pair at small 
momentum above the threshold in the corresponding partial wave. 

It can be noted that the constants C$ and Cp encode the information about the ampli- 
tudes for the bb pair to be in the corresponding color and orbital state as well as the amplitude 
for the conversion by the gluon operator E of the initial light quark-antiquark-gluon 'envi- 
ronment' into the final pion. Clearly, these constants are beyond present theoretical control, 
and for this reason it is not possible within the present approach to establish a quantitative 
relation between transitions to S- and P-wave states of bottomonium. The only guidance 
on the behavior of the light-hadron part of the amplitudes in Eq. ([2]) is provided by the soft- 
pion properties, which mandate the factor E n in the transitions to T(nS) and the factor p n 
in the transitions to hb(kP) |2J. Once these factors are accounted for as in Eq.(T5]), all the 
dependence on the excitation number of the specific final bottomonium state is contained in 
the overlap integrals ([3]) and (TjJ. 

In order to evaluate the latter overlap integrals we use the potential model of heavy 
quarkonium with the Cornell potential [9] 

V = ~- + - 2 (5) 

r cr 

with k = 0.52 and a = 2.34 GeV -1 , and calculate numerically the eigenfunctions R n s and 
RkP (we also set rrib = 5 GeV). We further consider the relative momentum q of the b and 
b quarks in the initial state as small. In the limit, where the Coulomb-like repulsion in the 
octet state is neglected, in the limit of small q the radial function in the 5*-wave state can 
be considered as constant -R^^r) const, while that in the P can be set as proportional to 
r: Rp{r) ~ const r. (Clearly, the overall normalization of these functions is not important 
for calculation of the ratios of the amplitudes A n s with different n and the ratios of A k p 
with k = 1 and 2.) In what follows we consider the modification of the overlap integrals in 
Eqs. ([3]) and (T4]) by the short- dsitance effect of the Coulomb-like repulsion between b and b 
in the color octet state. This repulsion is described by the potential 

Vs(r) = ^ , (6) 
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and we use in our estimates the value of the coefficient Kg related to that in the color singlet 
(Eq.flS])) as in a one gluon exchange: k 8 = k/8 = 0.065. 

One point related to the Coulomb-like repulsion that can be mentioned is that if the 
potential (E]) was applicable at all distances then the continuum wave functions would vanish 
in the limit q — > at any finite r due to impenetrability (from long distances) of the Coulomb 
barrier at zero energy. However neither the expression ([6]) is applicable at long distances, nor 
the momentum q is set to be literally zero. At longer distances the motion in the molecular 
state is described by that of heavy mesons, rather than individual heavy quarks, and also the 
typical values of q in the considered problem are small but finite and are set by the inverse size 
of the molecule. The overlap integrals are determined by the behavior of the wave function of 
the bb quark pair at short distances, i.e. at the characteristic size of bottomonium. At these 
distances the r dependence of the small- momentum wave function can still be calculated in 
the potential (EJ), while the normalization of the wave function is determined by the long- 
range modification of the interaction (j6]). Since the normalization of the functions cancels in 
the discussed here ratios of the amplitudes, one can use in a calculation of the integrals in 
Eqs. (EJ) and (jl]) either the small momentum limit of the continuum wave functions in terms 
of their short-distance expansion: 



(these terms do not depend on q), or introduce a small but finite q, and use the exact 
Coulomb functions (see e.g. in the textbook [TO] ) 



Rf = const e iqr X F X 1 + i-^] 2;-2iqr\ , = const r e iqr ^ 2 + i-^; 4; -2i q r 



with iFi(a;b;z) being the Kummer confluent hypergeometric function. We apply both ap- 
proaches and find that they result in similar estimates of the ratios of the considered overlap 
amplitudes. In particular we find that these ratios only weakly depend on q at q < 200 MeV, 
as is expected as long as q is not much larger than m& k 8 /2 w 160 MeV. 

The numerical results of our calculation and the experimental data are presented in the 
Tables 1 and 2. One can readily see that our estimates are within the range allowed by the 
current data. It is clear however, that there is much room for improvement of the data as 
well as for refinement of the theoretical approach. In particular, on the theoretical side, the 
specific numbers are fully dependent on the model wave functions for the bottomonium states 
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Table 1: Ratios of decay rates for Z&(10610) 



Ratio 


k 8 = 


Eqs. © 


Eqs. ( 


ED, q = 0^0.2 GeV 


Experiment [6] 


r 


Z b ^T(lS)n 




0.11 


0.09 




0.10 -=-0.11 


0.073 ± 0.029 


r 


Z b ^T(2S)n 






r[Z 6 ->T(35)7r 




0.62 


0.74 




0.70 -=-0.60 


0.49 ±0.19 


r[Z 6 ^T(25)7r 






T[Z b ^h b (2P)n] 
r{Z b ->h b (lP)n] 


0.58 


0.78 




0.72 -=-0.63 


1.54 ±0.95 


Table 2: Ratios of decay rates for Z{>(10650) 


Ratio 


K 8 = 


Eqs. © 


Eqs. ( 


I]), q = 4-0.2 GeV 


Experiment [6] 


r[Z£->T(is> 




0.10 


0.08 




0.09 4-0.10 


0.10 ±0.04 


r[Z^-s>T(25)vr 






r[zj->-T(35)7r 




0.86 


1.02 




0.97 4-0.83 


0.68 ±0.24 


r[Z^-+T(25)7r 






r[z^/i 6 (2P>j 
r[z^/i fc (ip>] 


0.73 


0.99 




0.91 4-0.80 


1.99 ± 1.11 



and on a general picture of the motion in a near threshold 'molecule'. Once more precise 
data might become available this may contribute to a better understanding of the structure 
of both bottomonium and of the molecular states of heavy mesons. Also, in our estimates we 
used a soft-pion approximation, and ignored any effects of a possible form factor depending 
on the momentum of the pion, which effects can be especially significant in the transitions 
to the final state T(lS)n. Such effects may arise from the unknown at present amplitude of 
the conversion of the light component of the meson-antimeson pair to pion (in the factors Cs 
and Cp) as well as from the recoil factors in the dipole matrix elements A n $ and A^p. The 
recoil factor is in fact determined by the process of conversion, namely by the fraction of the 
pion momentum transferred to individual heavy quark or antiquark as opposed to the recoil 
against the pair bb as a whole. (This is different from e.g. a photon emission, where the entire 
photon momentum is transferred to an individual quark or antiquark.) We are not aware at 
present of a proper way of including and estimating these momentum-dependent factors, and 
for this reason we chose to neglect these altogether. The fact that our numerical result for 
the relative yield of T(1S)tt is in a reasonable agreement with the data appears to indicate 
that the form factor effect should not be dramatic. One can also readily notice that our 
estimates for the yield of /i&(2P) relative to that of hb(lP) are about twice smaller than the 
central values of the experimental data and the agreement is only due to the currently large 
experimental uncertainty. If future more precise data would change this to a meaningful 
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disagreement, this could possibly point to a considerable deviation from the simple potential 
model used in our estimates, or indicate an enhanced contribution of dipole transitions to 
the P-wave bottomonium from a continuum .D-wave of the color-octet bb pair. 

MBV acknowledges numerous enlightening discussions of experimental and theoretical 
issues with Alexander Bondar. The work of MBV is supported, in part, by the DOE grant 
DE-FG02-94ER40823. 
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